Magnetic properties of the spin-| XXZ model on the Shastry-Sutherland lattice: 

Effect of long-range interactions 
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We study magnetic properties of the S = 1/2 Ising-like XXZ model on the Shastry-Sutherland 
lattices with long-range interactions, using the quantum Monte Carlo method. This model shows 
magnetization plateau phases at one-half and one-third of the saturation magnetization when ad- 
ditional couplings are considered. We investigate the finite temperature transition to one-half and 
one-third plateau phases. The obtained results suggest that the former case is of the first order and 
the latter case is of the second order. We also find that the system undergoes two successive tran- 
sitions with the 2D Ising model universality, although there is a single phase transition in the Ising 
limit case. Finally, we estimate the coupling ratio to explain the magnetization process observed in 
TmB 4 . 

PACS numbers: 75.40.Mg; 75.10.Jm; 75.40.Cx 
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Magnetic properties of quantum spin systems with 
frustrated antiferromagnetic interactions have been an 
interesting topic from both theoretical and experimen- 
tal aspects, because the strong fluctuations prevent the 
stabilization of classical orderings. Some exotic ordered 
states, such as a spin supersolid state on the triangular 
latticep], a spin liquid on the kagome latticey], 0] , 
and a t/(l) liquid state on the pyrochlorc lattice [5|, are 
good examples reflecting such fluctuations. The S — 
1/2 antiferromagnetic Heisenbergmodel on the Shastry- 
Sutherland lattices (SSLs)d, 0,11, has been studied 
as the other class with interesting characteristics derived 
from frustration and quantum fluctuation. In previous 
studies, it was suggested that this model has rich magne- 
tization phases including spin supersolid phases at mod- 
erate fillings between magnetization plateau phases [Toj. 
From experimental observations for the SSL compound 
SrCu(B03)2[8[, it was also discussed that there is the 
possibility of fractionalized magnetic plateaus at the 
magnetization m=m z /m s =l/<3, 1/6, 1/7, 1/9 and 2/9 in 
addition to the previously reported plateaus at m=l/3 
and l/4[9(. Here, m s indicates the saturation magnetiza- 
tion. While one theoretical scenario for explaining these 
plateaus was proposed in analogy to the quantum Hall 
effect @, a few points, such as the number of magnetiza- 
tion plateaus and the magnetization values of plateaus, 
are at the moment controversial. 

In recent experiments on rare-earth metals RB4 [R 
is a rare-earth atomjpj], [l2|, [l3|, 0, HH, fractionalized 
magnetic plateaus were also discovered at a very low 
temperature. These compounds have a tetragonal crys- 
tal structure P4/mbm and the magnetic moments orig- 
inating from R 3+ locate on the SSL in the afo-plane. 
In TmB4, a large m=l/2-magnetization-plateau region 
was confirmed for 1.9[T]< H <3.6[T] at a low temper- 
ature when the magnetic fields are applied along the 
c-axis[l3|, An important point that is the most 

different from SrCu(B03)2 is that, in the TmB4 case, 
a strong anisotropy along the c-axis is expected owing 
to the crystal fields. Recent analysies for specific heat 



measurements and the magnetization process [14j sug- 
gested that the J=6 multiplet of Tm 3+ is lifted, and 
then, the lowest energy state for a single ion is the 
non-Kramers doublet with J 2 =±6, having a large en- 
ergy gap between the lowest states. This leads us to 
describe the low-energy magnetic part of TU1B4 by a bi- 
nary model. Since it is expected that isotropic interac- 
tions work between the moments via itinerant electrons, 
we start from the two-site Hamiltonian, Hq + H', where 
Ho=-DJ2t(Jn 2 : H'=GJ2 l3 J* ■ Jj'j and 0<|G|<£>. In 
the ground state of H , there are four degenerated states, 
namely |J?,J?)=| ± 6, ±6} and | ± 6,=f6). When we 
treat H' as the perturbation, the first nontrivial off- 
diagonal term occurs from the 2 J-th order. The obtained 
effective Hamiltonian is the S — 1/2 Ising-like XXZ 
Hamiltonian with the transverse term proportionally to 
— {G/(— D)} 2J . Note that the transverse term always 
becomes ferromagnetic independently on the sign of G 
because J is integer. Meng and Wessel studied the mag- 
netization curves for the effective Hamiltonian by quan- 
tum Monte Carlo computations fill] . From the obtained 
results, it was found that the origin of the m=l/2 plateau 
can be explained by the quantum effect because only the 
m=l/3 plateau was confirmed in the Ising model case. 
However, the phase diagram for the Ising-like XXZ model 
case showed that the m=l/3 plateau phase spreads more 
widely than the m=l/2 plateau phase. Since the m=l/3 
plateau was not observed in TmB4, their results seem 
to be inconsistent with the experimental observations. 
This inconsistency may arise from the absence of "long- 
rang" interactions among magnetic moments of Tm 3+ 
ions, namely, the RKKY interactions. The magnetic 
properties of the Ising-like XXZ model with the long- 
range interactions have not studied yet. Thus, the clar- 
ification of the effects of long-range interactions on the 
magnetic properties is an important problem and helps 
us understand the origin of the m—1/2 plateau observed 
in T111B4. 

In this paper, we discuss the magnetic properties of 
the S = 1/2 Ising-like XXZ model on the SSL using 
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FIG. 1: (a) Effective model on the SSL with diagonal coupling 
7, nearest-neighbor coupling 7', and additional couplings J3 
and J4. The J3 coupling is defined on every plaquette with- 
out diagonal coupling, and the J4 coupling corresponds to 
the third-nearest neighbor coupling of the orthogonal-dimer- 
lattice geometry, (b), (c), and (d) are typical spin configu- 
rations in the m = collinear, f/3 plateau, and 1/2 plateau 
states, respectively. In (b)~(d), solid (open) circles represent 
up spins (down spins) and the broken (solid) line on each 
diagonal bond indicates an antiparallel (parallel) spin pair, 
(e) shows the bond Neel state in the intermediate phase (see 
text). 



the quantum Monte Carlo method based on the modi- 
fied directed- loop algorithm [17j. We consider the effects 
of J3 and J4 couplings in addition to the conventional 
SSL model with J and J' couplings. The geometrical 
configuration of the additional couplings J3 and J4 is 
shown in Fig. [T] (a) . The Hamiltonian considered here is 
described by 



w = E (JSi • Si) A± + E ( J ' s * • s ^)z 
+ e w • s 3 -) Ax + E • s,o, 

(*,i>" 
- SMs^E^ 2 ' 



(1) 



SiSj) + JS Z S Z 



where (JS< ■ S,-)^ = -A X |J| (S?S] 
and (ij), (ij)', (ij)", and (ij)"' mean sums over all pairs 
on the bonds with the J, J', J3, and J4 couplings, re- 
spectively. Note that the positive (negative) sign of each 
coupling means an antiferromagnetic (ferromagnetic) in- 
teraction and J > and J' > are considered. In the 
following computations, we set gjjiB = &b = 1 and treat 
the L x L system with the periodic boundary condition. 

We present the obtained magnetization processes in 
Fig. [5]when the easy-axis anisotropy and the coupling ra- 
tio are fixed at A^=0.2 and J' /J— 0.3, respectively. The 
magnetization-plateau states become stable at m = 1/2 
(to = 1/3) for 1.85< H/J <2.2 (1.25< H/J <1.55) 
when the weak ferromagnetic J4 (J3) coupling is taken 
into account. To investigate the spin configurations in 
both plateau states, we calculate two types of correla- 
tion functions. One is the bare spin correlation defined 
by C zz {r) = (S z (r)S z (0)) - (to z ) 2 , and the other is the 
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FIG. 2: Magnetization processes at ksT/J = 0.05. (a) and 
(b) are results when we consider J3 and J4 couplings, respec- 
tively. Circles (triangles) indicate the results for the addi- 
tional ferromagnetic (antiferromagnetic) couplings. The mag- 
netization processes without the additional couplings are rep- 
resented by squares. The symbols are drawn with error bars, 
which are smaller than the symbol size (here and the following 
figures) . 
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FIG. 3: Bare spin and bond spin correlations in 1/2 and 1/3 
plateau states at fcsT/7 = 0.05. The left-hand figures (al) 
and (bl) show the results for the 1/3 plateau state at J3/J = 
-0.106 and H/J = 0.75. The right-hand figures (a2) and (b2) 
show the results for the 1/2 plateau state at J4/7 = —0.106 
and H/J — 1. The horizontal axis is the distance along the 
(1,1) direction, which is normalized by the linear dimension 
L = 48. 



bond spin correlation defined by r^r^) = (B(r d )B(0)) — 
(B(r d )) 2 , where B(r d ) = S z rd+5 S z rd „ 5 and B(r d ) is de- 
fined only on the plaquettes with the diagonal coupling 
J. The bond spin correlation is convenient for character- 
izing the crystallization of the triplet dimers with S z — 1 
and S z = 0, which cannot be detected from the bare spin 
correlation. From Fig. [31 we find that both spin corre- 
lations show the presence of a true long-range order at 
to = 1/3 and m = 1/2. In the L x L lattice systems at 
UbT / J = 0.05, we confirm that no other plateau survives 
in the thermodynamic limit and that the ferromagnetic 
J4 coupling tends to stabilize only the 1/2 plateau even 
when the coupling ratio J'/J varies. Thus, we believe 
that the ferromagnetic J4 coupling plays a significant role 
in explaining the 1/2 plateau observed in T111B4. 

To discuss the reason why the field range of the 1/3 
or 1/2 plateau is expanded by the ferromagnetic J3 or 
J4 coupling, it is constructive to consider the Ising limit 
case. The lowest energy having the magnetization to = 0, 
1/3, and 1/2 can easily be estimated from each spin con- 
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figurations at T = 0. Fig. Q] shows the conventional spin 
configurations at a very low temperature, which are ob- 
tained from computations for the Ising model case. At 
m = 0, the collinear order shown in Fig. [1] (b) is stabi- 
lized by the J 3 and J 4 couplings, and the local energy is 
given by e co = — J/8 — J 3 /4 + J 4 /2. In the same manner, 
we obtain e 1/3 = -J/24 - J'/6 + J 3 /12 + J4/6 - H/6 
for the 1/3 plateau state and eu2 — J4/2 — /J/4 for the 
1/2 plateau state. Since the energy of the fully polarized 
state is given by e F = J/8 + J'/2 + J 3 /4 + J 4 /2-/J/2, the 
field ranges of the 1/3 and 1/2 plateaus can be calculated 
as A/J1/3 = 3 J' - 3 J 3 + 6 J 4 and AH 1/2 = 2 J 3 - 4 J 4 re- 
spectively. Consequently the ferromagnetic J4 (J 3 ) cou- 
pling makes the 1/2 (1/3) plateau state stable, which is 
in contrast to the antiferromagnetic case. 

Next, we investigate the finite temperature transition 
to the plateau phases. To assess the nature of the finite 
temperature transition to the 1/3 and 1/2 plateau states, 
we calculate the temperature dependence of the specific 
heat, and the results are shown in Fig. [U The spe- 
cific heat at (J 3 / J, J 4 /J, H/J) = (-0.106,0,0.725) has 
a single peak and the maximum value strongly diverges 
around T/J ~ 0.052 for increasing the system size L. 
This implies that the phase transition from the paramag- 
netic phase to the 1/3 plateau phase is of the first order. 
On the other hand, the weak system size dependence in 
the model at (J 3 /J, J 4 /J, H/J) = (0,-0.106,1.0) seems 
to result from the second-order phase transition with a 
negative exponent a < 0. 
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To clarify the order of the finite temperature transi- 
tion to the 1/3 plateau phase, we calculate static struc- 
ture factors both for bare spins S(Qo) and for bond 
spins r b (Q ) at Q = (7f/3,0). The results are shown 
in Fig. [5l Both structure factors suddenly develop at 
T/J ~ 0.052 and the values increase proportionally to 
the system size below the onset temperature. This re- 
sult also indicates the presence of the first-order tran- 
sition. In the phenomenological renormalization-group 
treatment for the first-order transition[l^], the bond-spin 
correlation ratio around the critical temperature scales 
as r b (J/2,L/2)/r b (J/4, J/4) - F(tL^ u ), where F is a 
scaling function, t = (T — T c )/T c , and the critical expo- 
nent v is given by the inverse of the space dimensionality 



l/d(d — 2). From Fig. [5] (c), we find that the finite-size 
scaling analysis is successfully performed and the data 
collapse is confirmed at T c /J = 0.516(3). 
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FIG. 5: (Color online) Transition to 1/3 plateau phase. Tem- 
perature dependence of static structure factors for (a) bare 
and (b) bond spins at Qo = (n/3,0). Insets show their size 
dependence, (c) Finite size scaling for correlation ratio of 
bond spins. 

We continue to study the universality class of the phase 
transition to the 1/2 plateau. From the low-temperature 
spin configuration in the Ising limit case (see Fig. Q](d)), 
we can expect C 4 symmetry breaking. Since the lowest 
energy state of the 1/2 plateau phase is fourfold degener- 
ated in context of the bare spin language, the universality 
class of this phase transition is naively expected to be the 
four-state-Potts universality class (the critical exponents 
are given by v — 2/3 and f3 — 1/12). In fact, we confirm 
a single phase transition belonging to the four-state Potts 
universality class in the Ising limit case[2l|. However, it 
is also possible that the C 4 symmetry breaks down in two 
separate steps as we see below. 

By introducing the 90 degree lattice rotation "c 4 " 
around the center of a plaquette without diagonal cou- 
pling, the symmetry group can be expressed as C 4 = 
{e, c 4 , c 4 , C4 1 }. At the higher critical temperature T c ±, 
the symmetry breaks down to {e, c|}, which is the sym- 
metry of the intermediate phase, and at the lower crit- 
ical temperature T C 2, the remaining symmetry breaks 
down to the trivial group {e}. To detect them, we in- 
troduce two kinds of order parameters here. The sym- 
metry breaking at T c \ can be characterized by a bond- 
spin "staggered magnetization" represented by B st = 
I ^ rd (— l)^( r<i ^i3(r ( i)|, where /(r<j) takes ±1 depending 
on the position of the diagonal coupling. The other sym- 
metry breaking at T C 2 can be characterized by freezing 
the antiferromagnetic fluctuation on the diagonal bond 
J, so that the order parameter is given by m c x = S% — S* 
or niy — S3 — S%. Here, the suffixes of longitudinal 
spin operators represent the site indexes shown in Fig. 
Q] (d). Fig. [6] (al) shows the temperature dependence 
of the Binder ratio Rb for the staggered magnetization 
B st . We find that Rb has a single crossing point at 
^bTci/J — 0.0822(3). Since the Z2 symmetry breaking 
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at T c x is expected from the viewpoint of bond spin order- 
ing, an intuitive speculation for the universality class of 
the phase transition leads us to expect the 2D Ising uni- 
versality class. We thus perform the data collapse assum- 
ing the scaling form T b (L/2, 0)/T b (L/3, 0) = F{tL 1 / v ) 
and BstL^I" = F(tL x l v ) with v = 1 and /3 = 1/8. 
The results are shown in Figs. [5] (a2) and (a3). The 
data collapse is confirmed with good accuracy despite 
the fact that we did not vary both the critical expo- 
nents and the critical temperature in the analysis. The 
temperature dependence of the Binder ratio Rs of m% 
provides a clear evidence of another phase transition. 
As shown in Fig. [6] (bl), Rs has a single crossing 
point at ksTc^l J = 0.0760(2) and the value is clearly 
different from that of the bond-spin ordering tempera- 
ture. The universality class of the phase transition at 
T C 2 is expected to be the same as that of the 2D Ising 
model because the remaining symmetry C2 is isomorphic 
to Z 2 - Actually, the finite-size scaling analysis for the 
fourth order cumulant U§ of m% can be achieved by us- 
ing fjj = F(tL 1 l v ) and v = 1, as shown in Fig. © (b2). 
From the obtained results, we conclude that there exists 
an intermediate phase that can be characterized by the 
bond Neel orderings accompanying the internal antifer- 
romagnetic bare-spin fluctuation (see Fig. Q] (e)). 

The obtained result has an interesting nature; the crit- 
icality is different from that of the four-state Potts model. 
One of the interpretations for this discrepancy is given 
by the phase diagram of the generalized four-state clock 
model, which may be similar to that of the generalized 
Ashkin- Teller modelpjj |2(j. In the Ising limit case, the 
phase transition corresponds to the P4 fixed point on the 
self-dual line. If the transverse coupling (5+ SJ~ + h.c.) 
is added, the criticality may change owing to the break- 
ing of self duality; the model undergoes two successive 
phase transitions with the 2D Ising universality or only 
one phase transition with the weak 2D Ising universal- 
ity. Further investigations in this direction are now in 
progress and will be published elsewhere (2l|. 

Finally, we discuss the magnetization of TmB4. It was 
observed that the 1/2 plateau state appears for 1.9[T]< 
H <3.6[T] at 2[K]43|, and the observed magnetization 
process associates with the coupling ratio J'/ J ~ 1 and 
the strong Ising anisotropy << l[l5[. From these 
facts, we roughly estimate the coupling ratios J3/J and 
J 4/ J by the local energy estimation in the Ising model 
case with J' / J = 1. We obtain the coupling ratios 
as J 3 /J ~ 0.1182 and J 4 /J < -0.25 at J'/ J = 1 
with J ~ 0.090 [K] [22| . Note that the field range of the 
1/2 plateau state does not depend on J 4 , although the 
J4 coupling is important for the appearance of the 1/2 
plateau. Furthermore we comment that this parameter 
set also gives a consistent Neel temperature ~ 7[K] at 
H = (experimentally, it was Tjy ~ 9.6[K]). In previous 
studies[la,[lg], it was clarified by the Monte Carlo calcu- 



lations for the minimal Ising model with J3 = J 4 = that 
the plateaus except the 1/3 plateau are not stabilized [l6j. 
Thus, our estimation seems to be reasonable. However, 
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FIG. 6: Transition to 1/2 plateau phase, (al) Temperature 
dependence of the Binder ratio Rb- Finite-size scaling analy- 
sis for (a2) correlation ratio and (a3) staggered magnetization 
B st . (bl) Temperature dependence of the Binder ratio Rs and 
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checking its validity more precisely via comparison with 
the other experimental observations might be required, 
because a small plateau at m ~ 1/8, which our model 
couldn't replicate it, has been observed experimentally. 
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